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Abstract 

Within the standard model there exist certain relations between 
CP violating rate differences in B decays in the SU(3) limit. We study 
SU(3) breaking corrections to these relations in the case of charmless, 
hadronic, two body B decays using the improved factorization model 
of Ref.@. We consider the cases B -» PP and B -» PV for both B d 
and B s mesons. We present an estimate for Acp(tt~tt + ) in terms of 
A C p(K-7t+). 

B decays are a subject of very active research at present since they provide 
useful information on the dynamics of strong and electroweak interactions for 
testing the Standard Model (SM) and models beyond and are ideally suited 
for a critical analysis of CP violation. The mixing induced CP asymmetry in 
B° —>■ ipK s versus B° — ► ipK s has already provided an accurate measurement 
of sin 2(3 [|l], 0. This result is in excellent agreement with the SM. Other mix- 
ing studies, such as B — ► 7r _ 7r + , are underway for determining sin 2a, but 
require more data to reduce the theoretical uncertainty associated with pen- 
guin contributions. Rate asymmetry measurements in the branching ratios of 
B mesons into mesons involving light quarks are also underway. These shed 
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light on direct CP violation in the amplitudes. Analyses of these decays to 
extract fundamental parameters of the SM are more difficult because of the- 
oretical uncertainty in the calculation of hadronic modes. In general, these 
asymmetries arise from interference of a Cabbibo suppressed tree amplitude 
with a (possibly enhanced) penguin amplitude. As such, these asymmetries 
are sensitive to contributions through loops, that could involve physics be- 
yond the SM. Thus, the study of direct CP violation can be a powerful tool 
to probe for physics beyond the SM if the theoretical uncertainty can be 
reduced. 

The goal of this letter is to study the direct CP violation asymmetry in 
a class of processes where there has been recent theoretical progress. These 
processes involve B decays into two light pseudoscalars P1P2 or into a light 
pseudoscalar and a light vector meson PV . We identify relations between 
rate asymmetries which are valid in the SU(3) limit in the standard model, 
and we compute SU(3) breaking corrections to them using the QCD improved 
factorization model of Ref . |§ . We also discuss additional relations which are 
valid in the SU(3) limit when annihilation contributions are neglected. 

CP violation in the SM arises solely from the phase in the 3x3 unitary 
CKM matrix, Vqkm = (Vij), and any CP violating observable is propor- 
tional to Im(VijV^V^Vki) , with i ^ k and j 7^ /. This simple property has 
important implications, as for example it leads to relations among CP vi- 
olating rate differences, Ap P = T(B — > PP) — T(B — > PP), for different 
decay modes. For instance, it has been shown that, with SU(3) flavor sym- 
metry, when small annihilation contributions and phase space differences are 
neglected, naive factorization yields the relation 

a 50 ~ \ B ° m 



This can be used to test the SM CP violation, or to predict one rate difference 
if the other one is known. The above equation leads to a relation for the CP 
violating rate asymmetry, 



Ac P ^)~-£j£^A CP (K-*+), (2) 



where Br(PP) are the CP averaged branching ratios, Br{ix 7r + ) = (5.2 ± 
0.6) x 1(T 6 and Br(K~-K + ) = (18.6 ± 1.1) x 1O _6 0, 0. Eq. g implies the fol- 
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lowing relation between the corresponding CP asymmetries: Acp{^ vr + ) pa 
-2AA CP (K-ir + ). 

Preliminary data on these asymmetries is just emerging, with BaBar re- 
porting A cp (tt-tt+) = -0.30 ±0.25 ±0.04, A cp (K-tt+) = -0.102 ± 0.050 ± 
0.0160 and Belle reporting A CP (n-n+) = 0.94±°;§f ± 0.09@. At the mo- 
ment the two experiments disagree on the value of Acp(7r~7r + ) but they still 
have very large errors. 

The most important question for theory is to establish the precision of 
Eqs. [I] and ^| within the standard model, or equivalently to estimate the 
corrections they receive. One can easily identify two sources of corrections 
for Eqs. [I] and 0: annihilation contributions, and SU(3)breaking effects. Even 
though the relation Eq. |I] already includes some SU(3) breaking effects in the 
factor f^/fx, it is necessary to have better control over these corrections in 
order to test the Standard Model. 

To begin our analysis of the B — > PP modes, we first note that there 
are several relations among the rate differences in these decays that follow 
from SU(3) flavor symmetry in the SM. There are also other relations such as 
Eq. [I] which rely both on SU(3) symmetry and on the neglect of annihilation 
amplitudes. It is easy to understand the origin of these relations. The decay 
amplitude for B — > PP can be parameterized as 

A{B -> PP) = V ub V: q T* P + V cb V c * q P^ (3) 

and can be decomposed into SU(3) invariant amplitudes according to the 
representation of the effective Hamiltonian H. SU(3) symmetry predicts 
that the amplitudes for B° — > 7r~7r + , K _ 7r + and B® — > tt~K + , K~K + are 
related and this can be proved by writing the decay amplitudes in terms of 
the SU(3) invariant amplitudes as 

T B° _ rpB° s - OAT , r T , r T , aT ,o r T 

Ti-« + = T!! K+ = C? + Cl -Al + 3C^, (4) 

where Ai indicate the annihilation contributions. Both model calculations 0, 
and fits to experimental data|| indicate that these annihilation amplitudes 
are small. The penguin amplitudes P PP can be parameterized in a similar 
way. 

We note that, even though T(P)f! 7r+ = T(P)^- K+ and T(P)f°_ w+ = 

B° 

T(P) n l K + in the SU(3) limit, there are no simple relations between the 
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branching ratios for these decays, because the CKM factors in the T and P 
amplitudes are different. However, because Ap P ~ Im(TP*)Im(V u bV* q V* b V cq ) 
and Im{V ub V* d V* b V cd ) = -Im(V ub V* s V; b V cs ) from the unitarity of the CKM 
matrix, we have the following relations among the CP violating rate differ- 
ences: 

A B° _ A m .BO _ A B° (t .s 

^7r-7r+ - ^K-K+i lX n-K+ ~ IX K^tt+ ■ 

These relations can be obtained by interchanging the d and s quarks (U-spin 
symmetry). If annihilation contributions are neglected, all the amplitudes 
T^°^ + , Tk-^ + , T^1 K+ and T^t K+ are approximately equal and one gets ad- 
ditional relations, 

A 7r-7r+ - -^K-K+ - /X tt-K+ - ~^K~-K+- y°) 

Similar relations exist as well for decays with neutral mesons in the final 
state. These relations are more complicated than those of Eq. |6] because 
there are more d and s quarks to interchange, and consequently it is harder 
to study the effect of SU(3) breaking in that case. For the remainder of this 
paper we concentrate on the relations in Eqs.[5] and ||. 

One must be careful, however, about the validity of these relations. In the 
exact SU(3) limit, the Ap P are also exactly zero because the standard model 
conserves CP when m s = ma- It is well known that in this case it is possible 
to remove the phase in the CKM matrix with an appropriate rotation among 
d and s quarks. In order to have a non-zero Ap P , one can not have an exact 
SU(3) symmetry. In order to have CP violation in the standard model no two 
quarks with the same charge can have the same mass; as long as m s ^ m& 
there will be CP violation regardless of how small these masses are. The 
relations in Eq. [| are thus valid and non-trivial in the limit where {m s — m^) 
is much smaller than the QCD scale, but not zero. 

When SU(3) breaking effects are included, the above mentioned relations 
will be modified and one needs a good understanding of these effects before 
using the relations to test the standard model. Our limited understanding of 
the strong interaction dynamics at low energies makes this task quite difficult. 
In what follows we illustrate the SU(3) breaking corrections that arise within 
the QCD improved factorization model of Ref. ||. 

Within this approach, the relevant decay amplitudes for B — > PP are 
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given byg, g] 



+ ^^(a^vrvr) + a? (vr7r) + rj(af (tttt) + ^(tttt)))] 

+ ^/B/.IWif") 

+ (K« + K^)fe(vrvr) + 26 4 (7T7r) 



S^{m\ - ml)F B -(m 2 K )f K [V ub V:MKn) 



+ V pb V p * s (al(Kir) + a? (Kvr) + r«(al(Kn) + al(Kn)))] 
+ i^f B Uf K [(V ub V: s + V cb V; s )(b 3 (K7r) - \bf w {Kv)l 

(7) 

where pis summed over u and c, r£ = 2m 2 /m b (m u +m ( i), = 2m 2 K j 'm b (m u + 
m s ) and 
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97TiV f 



^Mi(MM 2 ) + c 5 (A i 3 (MiM 2 ) + Al{M x M 2 )) 
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hiM^) = ^[c 4 A\(M 1 M 2 ) + c 6 Aj 2 (M 1 M 2 )}, 



bf w (M 1 M 2 ) 



6f w/ (M 1 M 2 ) 



= ( ^[c 9 A[(M 1 M 2 ) + c 7 (4(M 1 M 2 ) + 4(MxM 2 )) 

N c 

+ N c C s A f 3 (M 1 M 2 )}, 
C F 



[cwA\(M 1 M 2 ) + c 8 A 2 (M 1 M 2 )], 



(9) 



where Cf = (N% — l)/2N c and iV c = 3 is the number of colors and q are the 
Wilson coefficients. The vertex, the hard gluon exchange with the spectator, 
and the penguin contributions are: 



V, 



M 



pP 
r M,2 



pP 

r M,3 



pp,EW 
M,3 



12 In— -18+ f 1 dxg(x)§ M {x), 
a Jo 



Cl 



+ c 3 



H 

4 , rn b 2 . 
-ln- + --GM( Sp ) 

8 ln !^ + 4 Gm{0) _ Gm{1) 



+ (c 4 + c 6 ) 



^± i n _ ( n/ _ 2 )Gm(0) - G M (s c ) - G M (1) 



- 26 



eff 



1 



lo 1 — x 
P!iT = (ci + N c c 2 



$m(V), 



+ c 3 



3ln- + --G M ( Sp ) 



4 m 6 2 a / n 
3ln- + --G M ( Sp ) 

|ln — + 1-Gm(0)-G m (1) 



3c e 7 f r 1 -*e-* m ( X ), 

Jo 1 — X 



+ (c 4 + c 6 ) 



^ In ^ _ ( n/ _ 2)G M (0) - G M ( Sc ) - G M (1) 



- 26 



eff 



(ci + iV c c 2 ) 



4 77l 6 2 A / x 

3ln- + --G^( Sp ) 



— 3c' 



e// 
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H 



M\M2 



/b/tt 



1 dx 



+ 



m B X B F B ^(0) { Jo l-x $Ml(x) 



dy 
o l - y 



— ^Mi(x) 
X 



o l-y 



na s / dxdy FUx,y) , j = 1,3. 



A i J (M 1 M 2 ) 

A f 3 (M 1 M 2 ) = Kn n I (InhjEilv.u). 
with 

Ft(x,y) 

F^(x,y) 
Fl{x,y) 
Fi(x,y) 



(10) 
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(11) 



where x = 1 — x, y = 1 — y and the parameter 2{i M /m b coincides with r x . 
The functions g(x), Gm{%) and Gm(x) are given by 



g(x) = 3 
+ 

G(s,x) = 
G M (s) = 
G M (s) = 



l-2x 
1 — x 



In x — in 



2Li 2 (x) - ln 2 x + 

-4 



2 lnx 



1 — x 

— u) ln[s — u{\ — u)x] 

/ dxG(s — ze, 1 — x)$m(x), 
Jo 

/ dxG(s — ie, 1 — x)$ p (x 7 , 

JO 



— (3 + 2i7r) lnx — (x — > 1 — x) 



(12) 



where Sj = m\jm\ are the mass ratios for the quarks involved in the penguin 
diagrams, while <&m{x) and $p(x) are the distribution amplitudes of the M 
meson. The twist-3 distribution amplitude, $ p (x), is equal to 1, to the order 
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considered in the calculation. The distribution amplitude &m{x) has the 
following expansion in Gegenbauer polynomials ||. |7j 

$ M (x) = 6x(l - x)[l + a 1 C[ 3/2) {2x - 1) + a 2 C 2 3/2 (2x - 1) + ...], (13) 

with Ci(u) = 3u and C^ 2 (u) = (3/2)(5w 2 — 1), and is different for it and 
K. For 7r, the distribution in x must be even because the u and d quarks have 
negligible masses and their distributions inside the pion are symmetric. This 
dictates = 0. The coefficient a\ is estimated to be 0.1 ±0.3. For K, the u 
(or d) and s quarks inside the kaon are different, leading to an asymmetry in 
the x distribution. So a non-zero value for is needed and it is estimated 
to be 0.3 ± 0.3, while of = 0.1 ± 0.30, 0. 

One also has to consider divergences contained in the hard scattering and 
annihilation contributions. The divergent part in the hard scattering comes 
from Xh = Jq & p (x)dy / (1 — y) f=a Jq dy/(l — y), while the divergent part 
in the annihilation is of the same form at leading order. These divergences 
are logarithmic and, in principle, would be absent in a full theory. Here, we 
follow Ref. to introduce an infrared cut-off at A h = 0.5 GeV and use 

X H(A) =\n^. (14) 

The final results are insensitive to the precise value of the cut-off. As for 
the numerical inputs, we will use the values of the Wilson coefficients at 
fi = m 6 @, fi h = y/K^i, A h = 0.5 GeV, X B = 0.350 GeV. 

The decay amplitudes for B® — > tt~K + and — > K~ K + can be obtained 
from Eq. [7], by using the appropriate transition form factor F^ B ^ K and by 
changing l/m 2 B X B to l/m 2 B \ Bs in H Mi m 2 - 

Putting everything together, we are now able to estimate the size of dif- 
ferent contributions and, as expected, we find that the annihilation contri- 
butions are small. 

Eq. H] incorporates SU(3) breaking effects only through the difference in 
the decay constants f n and fx as they appear in naive factorization. To 
improve on this approximation we need to consider other sources of SU(3) 
breaking. For example there are SU(3) breaking mass differences in both the 
initial B mesons and in the final state mesons. These mass differences induce 
corrections that are proportional to m 2 M /m B and are therefore small. The 
decay amplitudes are proportional to the decay constant fu^ of Mi and to the 
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transition form factor F B ^ M2 , depending on which B is decaying into which 
final state. These form factors can also introduce SU(3) breaking effects. For 
B° — > 7r~7r + and B° — > K~tt + , all the corrections mentioned above account 
for the factor f%/ ' f\ in Eq. [|. Additional SU(3) breaking effects can arise from 
the distribution amplitude $>m(%), which is different for the d and s quarks. 
The important effect arises from the twist two distribution amplitudes. In 
our calculation we have used a constant $ p = 1 for the twist-3 distribution 
amplitude so we have neglected SU(3) breaking in this term. However, SU(3) 
breaking effects in this term should also be taken into consideration at higher 
order. 

To summarize, the SU(3) breaking effects that we do include are the 
difference in the decay constants and form factors; and the difference in the a\ 
and «2 terms that appear in the twist-2 distribution amplitude. With these 
effects taken into account, and using s u = = s s = and s c = (1.3/4.2) 2 in 
Eq. |H], the relation in Eq. |1] turns into 



7T~7T+ ^ J 7T 

/\B° ~ f2 



1 - 0.748a? - 0.109a£ - O.OUH^ - 0.004^ 
- 0.748«f - 0.109«f - 0.017^ + 0.00045f 



(15) 



where 5\ = 1-1.34X^-0.36(X5) 2 and 5$ = 0.1-0.8Xf +1.4(Xf ) 2 indicate 
the annihilation contributions. The numerical coefficients are obtained for 
the input parameters discussed before with X h ^a) real, and using X H ^) = 
h^ms/A). With these input parameters, H n7T and Hx n are in the range 
between 0.8 to 1. This leads to very small annihilation and hard scattering 



contributions as can be seen from Eq. 15. If one allows complex values for 
X H (A), then the corrections can be larger ||, but we do not have a good 
estimate for these parameters. 

The most important SU(3) breaking effect that we have identified (in 
addition to the difference in decay constants) arises from twist-2 distribution 
amplitudes. Using the central values a\ = 0, a\ = 0.1 and af- = 0.3, 
= 0.1, the size of A^_ 7r+ /A^._ 7r+ increases by a factor of 1.3, and is 
approximately —0.87. By taking into account the full range of values for the 
a-parameters, the maximum and minimum numerical values of the above 
ratios are, respectively, —1.4 and —0.6. This can be used to test the standard 
model and the QCD improved factorization. In particular, the sign of the rate 
difference is not changed by the SU(3) breaking effects we have considered. 

Similar calculations can be carried out for B s — > PP decays. We find 
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that the ratio of differences 



ABO ~ go 

is independent of the twist-2 distribution amplitudes or meson decay con- 
stants. This relation is particularly interesting because it can be used to 
test the SM with less uncertainties. The related CP asymmetries will be 
expressed in terms of the corresponding branching ratios which are scaled by 
transition form factors as 

Br(B» -> 7T-7T+) = C\i^\ Br(B»^n-K+)^, 

Ph~ 
~Ph 




Br(B°^K-n+) = C7 [ dbiH Br ( S s - K ~ K+ )^ , (17) 



KK 

where C = (m 2 B T B Jm 2 B T B ) and Phf iP2 = [(1- (m Pl +m ft ) 2 /m|)(l- (m Pl - 
rrip 2 ) 2 /m B } 1 / 2 /2m B is the phase space factor. In order to test the SM, one 
needs to know the form factors, these can be obtained from other processes 
or from theoretical calculations. Alternatively, one can use these relations to 
obtain the ratio of the form factors using experimental data. 

There are similar relations for B — > PV decays ||. By replacing one of 
the final octet pseudoscalar mesons with a corresponding octet vector meson 
in the previously discussed cases, one obtains 

^p-7r+ - K*~K+ ~ LA p-K+ — LA K—iv+ 

/\B° -_a^° -A 5 ? -_a«° nsi 

where the approximate sign indicates relations that hold true only when 
annihilation contributions are neglected. These relations, Eq. |TB|, are again 
expected to receive SU(3) breaking corrections. To estimate the SU(3) break- 
ing effects we use the QCD improved factorization model once again. 

When M 2 (the meson which picks up the spectator) is a vector meson, as 
for example in B° -> tt" p + , B° s -> K~K*+ and B° -> K~p+, B° s -> 7r-fT* + , 
the corresponding decay amplitudes can be obtained by replacing the form 
factor Fq^ p with Aq~* v and r x with — r x in Eq. [5], and by using the same 
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expressions for Eq. |9], except for Hm 1 m 2 which has no twist-3 terms. By 

is 



neglecting the annihilation contributions, the analogue of Eq. 15 



K- p + m B fl ( A^ p \ 1 + HOaf + 15.5af 



A B J. Rt+ rn Bs p K \A*^ K * ) 1 + llOaf + 15.5a? ' 



(19) 



and the same for A^L K * + /A^+. We observe the large coefficient of ai in 
both the numerator and denominator of Eq. IS. Since = and af = 



0.3 ± 0.3, the denominator has a vary large uncertainty, making a prediction 
for this asymmetry impossible within this framework. On the other hand, 
this provides an opportunity to constrain (or even to determine) af when 
the ratio in Eq. [19| is measured. 

When Mi is the vector meson, as in B° — > p~7r + , B° s — ► K*~K + and B® — > 
p~K + , B° — ■> K*~it + , the decay amplitudes can be obtained by replacing 

2 f ± 

the r x factor in Eq. [7] with r* K = w 0.3 (and similarly for r p ), and 

by removing the penguin terms Pj in the expressions for a 6 and a§ in 
Eq. H, because the vector meson is described only by a twist-2 distribution 
amplitude. With all this we obtain: 



Af-V _ m B / p 2 ( Ff~**\ 2 1 - 1.25a? - 0.18ag 



A S K L K+ rn Ba f*. \F^ K J 1 - 1.25af* - 0.18af * 



(20) 



Using the central values of the ranges a? = 0, a p 2 = 0.16 ± 0.09, af * = 
0.18 ± 0.05, af = 0.05 ± 0.05 and taking f p » 0.96/^0 we find: 

Our calculations show that important SU(3) breaking effects arise from 
the light-cone distributions of mesons in addition to those already present in 
the decay constants. These effects can only be estimated with large uncer- 
tainty because the parameters af 2 are not well determined at present. Using 
the currently allowed ranges we find, 

A CP (n-n + ) » - (3.1±5:|) A CP (K-n + ) , (22) 



1 We extract this ratio fp/fx ~ 0.96 from T(r — > p v)/T{t — » JT* v). 
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which can also be used to test the standard model and the improved factor- 
ization model to some extent. 

We have also shown that in the case of B — > PV, when the pseudoscalar 
meson is factored out, SU(3) breaking is large and estimates have very large 
uncertainty at present. In the case when the vector meson is factored out, 
as in Eq. the corrections are smaller. 

It is important to emphasize, however, that there are relations which 
are independent of a\ 2 parameters and decay constants. Examples include 
Eq. |l^, a corresponding relation for the ratio of branching ratios (Br(B° — > 
TT-TT + )/Br(B° -> K-ii + ) « Br(B s -> n-K + )/Br(B s -> K~K + )), and their 
analogues in B — > PV decays. These relations are more reliable than Eq. [I] 
in the sense that they do not receive the main SU(3) breaking corrections 
that we have investigated. Although this observation relies on the QCD 
improved factorization model, it may be more robust than model predictions 
for absolute values of rates because it only involves ratios. 

A systematic framework to study SU(3) breaking in B decays is, of course, 
needed before the relations we have presented can be used in precision tests 
of the standard model. With the estimates we have presented here, the 
relations are still useful. Large experimental violations of them would hint 
at possible new physics; at the very least they would provide information on 
the limitations of the QCD improved factorization model. To test some of 
the relations that we have discussed, charmless hadronic two body B s decays 
must also be measured. This can not be done by the B-factories at present, 
but in the near future such B s decays will be studied at the Tevatron II and 
at LHCb. 
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